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A  SIMPLIFIED  METHOD  OF  ELASTIC-STABILITY 
ANALYSIS  FOR  THIN  CYLINDRICAL  SHELLS 
II  -  MODIFIED  EQUILIBRIUM  EQUATION 
By  S.  B.  Batdorf 


SUMMARY 


A  modified  form  of  Donnell's  equation  for  the  equilibrium  of 
thin  cylindrical  shells  is  derived  which  is  equivalent  to  Donnell's 
equation  hut  has  certain  advantages  in  physical  interpretation  and 
in  ease  of  solution-,  particularly  in  the  case  of  shells  having 
clamped  edges.  The  solution  of  this  modified  equation  by  means  of 
trigonometric  series  and  its  application  to  a  number  of  problems 
concerned  with  the  shear  buclcling  stresses  of  cylindrical  shells 
are  discussed.  The  question  of  implicit  boundary  conditions  also 
is  considered. 


INTRODUCTION 


During  a  general  theoretical  investigation  of  the  stability 
of  curved  sheet  under  load,  a  method  of  analysis  was  developed 
which  appears  to  be  simpler  to  apply  than  those  in  general  use. 

The  development  of  this  method  is  presented  in  two  parts,  of  which 
reference  1  is  the  first  and  the  present  paper  the  second.  The 
specific  problems  solved  by  this  new  method  are  treated  in  detail  in 
other  papers.  (See,  foh  example,  references  2  to  7.) 

In  reference  1.  the  stability  of  a  stressed  cylindrical  shell 
was  analyzed  in  terms  of  Donnell's  equation,  a  partial  differential 
equation. for  the  radial  displacement  w,  which  takes  into  account 
the  effects  of  the  axial  displacement  u  and  the  circumferential 
displacement  v.  Reference  1  shows  the  manner  in  which  this  equation 
can  be  used  to  obtain  relatively  easy  solutions  to  a  number  of 
problems  concerning  the  stability  of  cylindrical  shells  with  simply 
supported  edges.  The  results  of  the  solution  of  this  equation  were 
shown  to  take  on  a  simple  form  by  the  use  of  the  parameter  k 
(similar  to  the  buckling— stress  coefficients  for  flat  plates)  to 
represent  the  state  of  stress  in  the  shell  and  the  parameter  Z 
to  represent  the  dimensions  of  the  shell,  where  Z  is  defined  by 
the  following  equations: 
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For  a  cylinder  of  length  L 

Z  =  l£  j  1  -  u2 
rt  ' 

and  for  a  curved  panel  of  width  b 


where 

r  radius  of  curvature 
t  thickness  of  shell 
and 

(a  Poisson’s  ratio  for  material 

The  accuracy  of  Donnell’s  equation  was  established  by  comparisons 
of  the  results  found  by  its  use  with  the  results  found  by  other 
methods  and  by  experiment. 

In  the  simplest  method  that  has  been  found  for  solving 
Donnell’s  equation,  the  radial  displacement  w  is  represented 
by  a  trigonometric  series  expansion.  This  method  can  be  used  to 
great  advantage  for  cylinders  or  curved  panels  with  simply  supported 
edges  but  leads  to  incorrect  results  when  applied  uncritically 
to  cylinders  or  panels  with  clamped  edges. 

In  the  present  paper  an  equation  is  derived  which  is  equivalent 
to  Donnell’s  equation  but  is  adapted  to  solution  for  clamped  as 
well  as  simply  supported  edges  by  means  of  trigonometric  series. 

This  modified  equation  retains  the  advantages  of  Donnell’s  equation 
in  ease  of  solution  and  simplicity  of  results.  The  solution  of  the 
modified  equation  by  means  of  the  Galerkin  method  is  explained,  and 
the  results  obtained  by  this  approach  in  a  number  of  problems 
concerned  with  the  shear  buckling  stresses  of  cylindrical  shells 
are  given  in  graphical  form  and  discussed  briefly.  Boundary  conditions 
implied  by  the  method  of  solution  of  the  modified  equation  are 
also  discussed. 
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SYMBOLS 


a 

L 

a±:>  *u 

ke 


P 

r 

t 

u 


w 


length  of  curved  panel  (longer  dimension) 
width  of  curved  panel  (shorter  dimension) 
deflection  coefficients  in  trigonometric  series 


shear— stress  coefficient 


TtV 

Drf2 


for  curved  panel 


for  cylinder; 


fax tLg 


direct-axial -stress  coefficient  { ——  for  cylinder; 

Vr^2 


OS rth2 

-  -■■■■  for  curved  panel! 

Drt2  y 


circumferential— stress  coefficient 

Cyth2  '  '  \ 

~'7  — g—  for  curved  panel  ! 


rr 

a—  for  cylinder 
^  Djt2 


Dit£ 

lateral  pressure 

radius  of  curvature  of  cylindrical  shell 

thickness  of  cylindrical  shell 

displacement  in  axial  (x— )  direction  of  point  on  shell 
median  surface 

displacement  in  circumferential  (y— )  direction  of  point 
on  shell  median  surface 

displacement  in  radial  direction  of  point  on  shell 
median  surface;  positive  outward 

axial  coordinate 

circumferential  coordinate 


n,  P,  <Lj 
D 


integers 

plate  flexural  stiffness  per  unit  length 


/  Et3  \ 

\I2(1  -  V2)) 
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E 

F 


Young’s  modulus  of  elasticity 


L 

Q»  Qi>  Q2 
Rs 

Ex 


P 

Ox 

CY 

T 

cr 


by  buckle  deformation 
d2F 


Airy’s  stress  function  for  median— surface  stresses  produced 

(  — =■,  stress  in  axial  direction; 

Vdy* 

S^F 

•,  stress  in  circumferential  direction;  -  .  shear 

ax2  Sx  Sy 

stress  } 

length  of  cylinder 
mathematical  operators 

shear— stress  ratio;  ratio  of  shear  stress  present  to 
critical  shear  stress  when  no  other  stress  is  acting 

axial-compressive-stress  ratio;  ratio  of  direct  axial 
stress  present  to  critical  compressive  stress  when 
no  other  stress  is  acting 

curvature  parameter  j/l  -  jjl2  for  cylinder; 

p  v 

~  j/l  —  |j2  for  curved  panel  or  long  curved  strip) 

half  wave  length  of  buckles;  measured  circumferentially 
in  cylinders  and  axially  in  long  curved  strips 

Poisson’s  ratio  for  material 

applied  axial  stress.,  positive  for  compression 

applied  circumferential  stress,  positive  for  compression 

applied  shear  stress 

critical  shear  stress 


operator 


a"  +  ^  a*  .  s4 


‘S?+a£^  +  i X 


r 
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v8 

V1'- 


inverse  operator  defined  by  equations 
V^O^f)  =  ^(V-^)  =  f 


THEORY 

Derivation  of  Modified  Equation 


The  equation  of  equilibrium  for  a  flat  plate  may  be  written 


IfA* 


+  tfcr 


8^w 

dyfi 


+  2r 


a  &V  p  = 

8x  by  y  by2/ 


(U 


where  p  is  lateral  pressure.  (This  equation  is  equivalent  to 
equation  (197)  of  reference  8.) 

For  a  cylindrically  curved  plate  having  a  radius  of  curvature  r, 
the  following  pair  of  simultaneous  equations  of  equilibrium  may  be 
written  (as  a  generalization  of  equations  (11)  and  (10)  of  reference  9) : 


DV,+w 


+  2t 


b^w 
bx  by 


t  8%  _ 
bx^ 


r 


(2) 


9^f  —  S  8^w  _  0 
r  8x2 


(3) 


where  F  is  Airyrs  stress  function  for  the  median— surface  stresses 
produced  by  the  buckle  deformation  (reference  10) .  Equation  (2) 

differs  from  equation  (1)  only  in  the  addition  of  the  term  —  , 

r  bx2  ' 

which  expresses  the_effect  of  the  curvature.  Equation  (3)  shows 
that.,  unlike  flat  plates,  cylindrical  shells  experience  stretching 
of  the  median  surface  when  originally  straight  lines  in  the  surface 


V 
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are  "bent  slightly.  Elimination  of  F  between  equations  (2) 
and  (3)  by  suitable  differentiations  and  additions  gives  the 
following  single  equation  in  w  for  the  equilibrium  of  cylindrical 
shells : 


DV8v  + 


5gw 

dx^ 


+  2t 


a2- 


w 


dx  dy 


+  av 


4 

+  Vj 


(4) 


Equation  (4),  which  was  first  derived  by  Donnell  (reference  11), 
was  treated  in  reference  1. 


An  alternative  method  for  obtaining  a  single  equation  in  w 
for  the  equilibrium  of  a  cylindrical  shell  is  to  solve  equation  (3) 
for  F  and  substitute  the  result  into  equation  (2).  This  procedure 
can  readily  be  carried  out  in  the  following  manner.  Differentiation 
of  equation  (3)  twice  with  respect  to  x  gives 


V1*  &§£  _  2  bjf w 

r  5? 


0 


The  symbolic  solution  of  equation  (5)  for 


dgF  E  \y  -4 


(5) 


Substitution  of  this  result  into  equation  (2)  gives 


Equation  (6)  is  simply  equation  (4)  modified  by  multiplication  by 
the  operator  v~4.  In  the  present  paper,  equations  (4)  and  (6) 
are  referred  to  as  Donnell’s  equation  and  the  modified  equation, 
respectively . 


Advantages  of  Modified  Equation 

One  of  the  quickest  and  most  convenient  methods'  for  obtaining 
solutions  of  flat— plate  buckling  problems  to  any  desired  degree  of 
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approximation  uses  a  Fourier  series  type  of  expansion  for  the  deflection 
surfaoe  w.  Both  Donnell  'a  equation  end  the  modified  equation  can 
he  solved  hy  this  method  in  the  case  of  buckling  problems  involving 
curved  plates  having  simply  supported  edges. 

As  mentioned  in  the  ‘'Introduction,"  however,  Donnell's  equation 
is  not  well  adapted  to  solution  by  Fourier  series  of  problems 
involving  the  stability  of  shells  with  clamped  edges .  The  cause  of 
the  trouble  appears  to  be  that  the  calculation  of  some  of  the  high- 
order  derivatives  found  in  Donnell 's  equation  Borne times  leads  to 
divergent  trigonometric  series  when  the  edges  are  clamped.  The 
modified  equation,  however,  is  applicable  to  clanrped-edge  problems 
as  well  as  to  problems  involving  simply  supported  edges  because 
lower-order  derivatives  are  involved. 

In  addition  to  its  advantages  in  the  solution  of  problems 
involving  shells  with  clamped  edges,  equation  (6)  has  the  additional 
advantage  that  each  term,  has  a  definite  physical  significance: 

The  first  term  gives  the  restoring  force  per  unit  area  of  the  deflected 
surface  due  to  bending  and  twisting  stiffnesses;  the  second  term 
gives  the  restoring  force  per  unit  area  due  to  Btretching  stiffness; 
and  the  remaining  terms  give  the  deflecting  forces  per  unit  area 
due  to  applied  loads.  Because  of  these  advantages,  the  modified 
equation  waB  adopted  for  general  use  in  references  2  to  7 . 

Both  Donnell's  equation  and  the  modified  equation  result  in 
the  Bame  critical  stresses  for  simply  supported  cylindrical  shells, 
and  the  two  methods  require  essentially  equivalent  mathematical 
processes.  (See  appendix.)  The  characteristics  of  solutions  by 
means  of  Donnell  'b  equation  in  the  case  of  simply  supported  BhellB 
(reference  1)  -  namely,  the  theoretical  cylinder  parameters,  the 
simplicity  of  calculations  and  results,  and  the  implied  boundary 
conditions  on  u  and  v  -  are  characteristics,  also,  of  solutions 
by  means  of  the  modified  equation.  The  same  characteristics,  except 
for  a  change  in  the  implied  boundary  conditions  on  u  and  v,  also 
apply  to  solutions  of  clamped-edge  shell  problems  by  means  of  the 
modified  equation.  ThiB  change  is  discussed  in  the  section  entitled 
'Boundary  Conditions . " 


Solution  of  Modified  Equation  by  Galerkin  Method 

An  approximate  method  of  solving  vibration  and  buckling  problems 
closely  paralleling  that  of  Bitz  was  introduced  in  1915  by  Galerkin. 
(See,  for  example,  references  12  and  13  •)  The  main  distinction  between 
the  Bitz  and  Galerkin  methods  is  that  the  Bitz  method  begins  with  an 


8 


NACA  IN  No .  1342 


energy  expression j  whereas  the  Galerkin  method  begins  with  an 
equation  of  equilibrium.  The  Galerkin  method  is  readily  adaptable 
to  the  solution  of  equation  (6)  and  is  now  desoribed  briefly. 

Let  the  equation  of  equilibrium  be  written 


Q(w)  +  p  =  0  (T) 

where  p  is  lateral  pressure  and  Q  is  some  operator  in  x  and  y 
which  for  the  purposes  of  this  paper  is  taken  to  be  linear.  According 
to  the  Galerkin  method.,  the  equation  may  be  solved  by  expanding 
the  unknown  function  w  in  terms  of  a  suitable  set  of  functions 
fi(x)gj(y),  each  of  which  satisfies  the  boundary  conditions  but 
not  in  general  the  equation  of  equilibrium: 


w 


=  T5 


/  ai,Jfi  00gj(y) 
i  J 


(9) 


Substitution  of  this  expression  for  w  into  equation  (7)  gives 
the  following  equation: 


ais5Qffi(x)g()(y)]  +  p  =  0  (9) 

i  J  ' 


Because  the  functions  fi(x)gj(y)  were  chosen  to  satisfy  the 

boundary  conditions  rather  than  the  equation  of  equilibrium, 
equation  (9)  cannot,  in  general,  be  satisfied  Identically  by  any 
choice  of  the  coefficients  &i.j.  These  coefficients  can  be  chosen, 

however,  to  assure  the  vanishing  of  certain  weighted  averages  of 
the  left-hand  side  of  equation  (9) .  The  weighting  functions  used 
in  the  Galerkin  method  are  the  original  expansion  functions,  so 
that  the  following  simultaneous  equations  for  determining  the 
coefficients  ap  ^  are  obtained: 

2_  ^9mniJaiJ  =  0  (10) 

i  5 

(m=l,2, 3,  •  •  •  >  n=l,2, 3,  •  •  • ) 


r 
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vhers 


■^mni,) 


fm(*>8n(y) 


<^Q[?i(x)gj(y)_ 


dx  dy 


(ID 


The  simultaneous  set  of  linear  algebraic  equations  in  the  unknown 
coefficients  ajj  (equation  (10)),  obtained  by  using  the  original 

expansion  functions  as  weighting  functions,  is  ordinarily  the  same 
set  which  would  be  found  by  the  Eitz  method,  if  the  same  series 
expansion  for  w  were  used.  A  solution  of  any  desired  degree  of 
accuracy  may  therefore  be  obtained  by  the  Galerkin  method. 

In  applying  the  Galerkin  method  to  equation  (6)  by  use  of 
Fourier  series  expansion  for  w,  expressions  of  the  type 


must  be  evaluated.  The  operator  V  the  inverse  of  V  , 
simply  introduces  into  the  denominator  of  each  tepn  of  the  series 
the  expression  that  comes  into  the  numerator  if  7’  is  applied. 
Thus, 


sin  (12) 

D 


This  result  may  readily  be  verified  by  applying  the  operator  V  ^ 
to  each  side  of  equation  (12) . 

In  writing  equation  (12)  the  quantity  V“^f,  as  defined  by 
the  equation 
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was  tacitly  assumed  to  be  unique.  The  quantity  actually  is  not 
unique ;  any  number  of  terms  which  vanish  when  operated  upon  by  V  ^ 
could  be  added  to  the  right-hand  side  of  equation  (12) .  The 
omission  of  such  terms  makes  the  present  analysis  parallel  to  the 
analysis  using  Donnell’s  equation  (see  reference  1)  and  implies 
certain  boundary  conditions  on  u  and  v,  which  are  discussed 
in  a  subsequent  section  entitled  "Boundary  Conditions.", 


Deflection  Functions 

Simply  supported  edges.-  For  simply  supported  cylindrical 
shells,  the  following  series  expansions  for  w  may  be  used  to 
represent  the  buckle  deformation  to  any  desired  degree  of  accuracy 
(in  these  functions,  x  is  the  coordinate  in  the  axial  direction 
and  y,  the  coordinate  in  the  circumferential  direction) : 

(1)  Rectangular  curved  plate  (axial  dimension  a  and 
circumferential  dimension  b) 

00  CO 

\. 

w  =  ^  y  sin  sin  (13) 

m-1  n-1 


(2)  Curved  strip  long  in  the  axial  direction  (circumferential 
width  b  and  axial  wave  length  21) 

(a)  Direct  stresses  only 


w  =  sin  23  \  a  gin  MI 

A.  Z _  b 

m=l 


(14) 


(b)  Shear  stress  with  or  without  addition  of  direct 
stress 


°°  00 

w  =  sin  £3^  am  sin  22Z  +  c08  Z& fcm  sin  (15) 

^1  fcr 


r 
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(3)  Complete  cylinder  (length  L  and  circumferential  wave 
length  2X) 

(a)  Direct  stresses  only 
00 

;Y 


w  =  sin  )  am  sin 


m=l 


mrtx 

L 


(16) 


(h)  Shear  stress  with  or  without  addition  of  direct 
stress 


CO 


(17) 


Clamped  edges.—  Probably  the  simplest  method  of  treating 
cylindrical  shells  with  clamped  edges  is  to  employ  the  expansions 
in  equations  (13)  to  (17)  modified  by  substituting  functions  of  the 
type 


<Pm(x)  =  sin  sin 


cos 


(m  —  l)nx 


cos 


(m  +  l)nx 
a 


(18) 


wherever  functions  of  the  type  sin  appear,  with  a  similar 

Q 

substitution  for  functions  of  y  (all  terms  involving  summation 
subscripts  m  and  n  are  thus  changed;  terms  involving  X, 

such  as  sin.  — ,  remain  unchanged) ,  The  functions  qxri(3c)  form 

a  complete  set  so  thatfinite  expansions  for  w  of  the  type  suggested 
for  shells  with  clamped  edges  as  well  as  those  for  shells  with  simply 
supported  edges  may  be  used  to  represent  the  buckle  deformation 
to  any  desired  degree  of  accuracy. 


Boundary  Conditions 

Simply  supported  edges.—  Appendix  D  of  reference  1  shows  that, 
if  the  buckling  stress  of  a  simply  supported  shell  is  found  by  means 
of  the  expansions  for  w  given  in  the  preceding  section  entitled 
"Deflection  Function,"  the  boundary  or  edge  conditions  implied  for 
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the  median-surface  displacements  u  and  v  are  zero  displacement 
along  each  of  the  edges  of  a  cylinder  or  curved  panel  and  free 
displacement  normal  to  each  edge .  (Although  the  proof  given  used 
equation  (4),  the  proof  could  equally  well  have  teen  based  on 
equation  (6).) 


The  boundary  conditions  for  simple  support  may  thus  be  -written, 
at  a  curved  edge  (x  =  constant) , 


^  5% 

5x2  5y2 


0 


(19) 


and,  at  a  straight  edge  (y  =  constant). 


w 


5^w  5^F 

5y2 


0 


(20) 


Clamped  edges.  -  By  a  method  similar  to  that  in  appendix  D 
of  reference  1,  solutions  using  the  functions  suggested  in  the 
preceding  section  for  the  treatment  of  clamped  edges  can  be  shown 
to  correspond  to  the  boundary  conditions  -  zero  displacement  normal 
to  an  edge  and  free  displacement  along  an  edge  • 

The  boundary  conditions  for  clamped  edges  therefore  become, 
at  a  curved  edge  (x  =  constant). 


w 


5w 

s 


=  u  = 


5x2 


and,  at  a  straight  edge  (y 


=  constant). 


w 


5 w  5fP 

3y  5y2 


0 


(22) 


Discussion. -  As  mentioned  in  reference  1,  the  boundary  conditions 
implied  for  u  and  v  in  the  case  of  simply  supported  edges  are 
appropriate  for  cylinders  or  panels  bounded  by  light  bulkheads  or 
deep  stiffeners,  which  are  stiff  in  their  own  planes  but  may  be 
readily  warped. out  of  their  planes. 
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The  "boundary  -conditions  on  u  and  v  appropriate  for  a 
clamped  edge  would  seem  to  "be  zero  displacement  normal  to  the  edge 
and  zero,  rather  than  free,  displacement  along  the  edge.  Comparison 
of  critical  stresses  for  shells  with  clamped  edges  found  "by  the 
method  in  the  present  paper  with  critical  stresses  found  "by  the 
method  in  references  9  and  1^,  giving  "boundary  conditions  u  a  v  *»  0, 
however,  indicates  that  the  imposition  of  the  added  requirement 
of  zero  displacement  along,  the  edge  ordinarily  has  very  little 
effect  on  the  critical  stresses. 

A  less  satisfactory  method  for  solving  problems  concerning 
shells  with  clamped  edges  involves  the  use  of  functions  of  the  type 


1 

m 


ein  _ 1-* 

a  m  +^2 


sin 


(m  +  2)itx 
a 


instead  of  those  described  by  equation  (l8) .  In  this  method,  the 
functions  used  are  those  for  simple  support  taken  in  such  combinations 
that  the  edge  slope  is  zero.  Use  of  such  functions  leads  to  the 
same  boundary  conditions  on  u  and  v  as  were  described  for  simply 
supported  edges;  at  the  edge  y  =  constant,  for  example,  the  boundary 
conditions  become 


w 


dgF 

8x2 


0 


(23) 


The  use  of  these  functions  to  represent  shells  with  clamped  edges 
is  not  recommended,  however,  for  the  following  reasons:  The  associated 
boundary  conditions  seem  to  be  artificial  and  unlikely  to  "be  reproduced 
even  approximately  in  actual  construction;  the  method  leads  in  some 
oases  to  solutions  that  differ  considerably  from  the  solution  for 
ideal  clamped -edge  conditions  in  which  u  =  v  =  0;  and  the  solutions 
obtained  generally  converge  rather  poorly. 


APPLICATIONS  AND  DISCUSSION 


Among  the  more  difficult  she 11 -stability  problems  to  treat 
theoretically  are  those  which  Involve  shear  stresses.  In  fact, 
until  I93I+,  when  Donnell’s  paper  on  critical  shear  stress  of  a 
cylinder  in  torsion  was  published  (reference  11),  such  problems 
wears  generally  regarded  as  impracticable  to  solve .  In  order  to 
Illustrate  the  type  of  solution  to  be  found  by  the  method  of  analysis 
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Just  outlined  and  the  effect  of  boundary  conditions  on  critical 
stresses,  the  results  obtained  for  a  number  of  shell  stability 
problems  involving  shear  stresses  are  reproduced  and  discussed 
briefly  here.  The  problems  treated  are  summarized  in  table  I. 

Critical  shear  stress  of  long  curved  strip. -  The  critical 
shear  stress  for  a  long  plate  with  transverse  curvature  is  given 
by  the  equation 


T 


cr 


k 


rt%> 

®b2t 


where  3%  is  a  dimensionless  coefficient,  the  value  of  which 
depends  upon  the  dimensions  of  the  strip,  Poisson 's  ratio  for 
the  material,  and  the  type  of  edge  support.  In  figure  1  (fig.  1 
of  reference  2)  the  shear -stress  coefficient  ks  is  given  for 
plates  with  simply  supported  edges  and  with  clamped  edges .  This 
solution  for  simply  supported  edges  coincides  with  that  given  by 
Kromm  (reference  15). 

As  indicated  in  the  previous  section  entitled  'Poundary 
Conditions,  "  the  solution  corresponding  to  the  boundary  conditions 
of  equation  (23)  (dashed  curve  of  fig.  1)  is  poorly  convergent  and 
deviates  appreciably  frcm  the  results  for  completely  fixed  edges. 

Figure  1  shows  this  poor  convergence  in  the  limiting  case  of  a  flat 
plate,  for  which  the  critical  stress  is  independent  of  boundary 
conditions  on  u  and  v.  Even  a  tenth -order  determinant  led  to 
a  result  that  is  7  percent  above  the  true  solution;  whereas  the 
result  using  a  fourth -order  determinant  obtained  with  the  deflection 
functions  recommended  for  clamped  edges  is  only  1  percent  above. 

In  figure  2  (fig.  2  of  reference  2)  the  solutions  given  in 
figure  1  are  compared  with  the  results  given  by  Leggett  (reference  9) 
for  simply  supported  and  clamped  edges  with  u  =  v  =  0  at  each 
edge.  Throughout  the  range  for  which  they  are  given,  Leggett's 
results  for  clamped  edges  differ  only  slightly  from  those  of  the 
present  paper.  On  the  other  hand,  the  previously  mentioned 
discrepancy  between  the  results  for  completely  fixed  edges 
(u  =  v  =  0)  and  those  for  the  boundary  conditions  of  equation  (23) 
(dashed  curve)  may  be  inferred  from  this  figure  to  be  considerable 
for  large  values  of  Z.  A  minimum  measure  of  this  discrepancy  is 
the  distance  between  the  clamped -edge  curves  for  v  =  0  and  for 
u  =  0  in  figure  2,  since  Leggett's  curve  must  alvayslie  above  the 
curve  for  v  =  0. 

The  reason  for  the  marked  increase  in  buckling  stress  of  simply 
supported  curved  strips  when  the  edges  are  restrained  against 
circumferential  displacement  during  buckling  is  discussed  in  reference  2. 
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Critical  shear  stress  of  cylinder  in  torsion. -  The  critical, 
shear  stress  of  a  cylinder  subjected  -to  torsion  is  given  by  the 
equation 


r 


cr 


=  k, 


Vt 


In  figure  3  (fig-  1  of  reference  3)  the  values  of  ks  are  given 
for  cylinders  with  simply  supported  edges  (boundary  conditions 
of  equation  (19))  and  cylinders  with  clamped  edges  (boundary 
conditions  of  equation  (21)).  At  high  values  of  Z ,  the  values 
of  ks  for  thick  cylinders  are  given  by  special  curves  for  various 

values  of  ^  \j  1  -  as  discussed  in  reference  1.  At  values 

of  Z  greater  than  about  100  only  a  small  increase  in  buckling 
stress  is  caused  by  clamping  the  edges .  The  results  indicated 
in  figure  3  are  in  very  close  agreement  with  Donnell's  results  for 
the  same  problem,  except  in  the  range  5  <  Z  <  500  vheye  the 
somewhat  lower  curves  of  the  present  paper  represent  a  more  accurate 
solution. 

Reference  1  shows  that  boundary  conditions  imposed  upon  u 
and  v  at  the  curved  edges  of  a  panel  or  cylinder  have  an  almost 
insignificant  effect  on  the  buckling  stresses,  whereas  conditions 
imposed  on  v  at  the  straight  edges  may  be  quite  important. 
Comparison  of  figure  1,  in  which  boundary  conditions  on  straight 
edges  are  considered,  with  figure  3 >  in  which  conditions  on  curved 
edges  are  considered,  indicates  that  a  similar  situation  exists 
with  respect  to  restraint  against  edge  rotation. 

Critical  shear  stress  of  curved  panel. ~  The  values  of  ke 
giving  the  critical  shear  stresses  of  simply  supported  curved 
rectangular  panels  are  given  in  figures  4  and  5  (figs*  1  nnd  2, 
respectively,  of  reference  4) •  The  corresponding  boundary  conditions 
on  u  and  v  are  zero  displacement  parallel  to  the  edges  and  free 
warping  normal  to  the  edges •  Figure  4  indicates  that  as  the 
curvature  parameter  Z  increases,  the  critical  shear  stresses  of 
panels  having  a  circumferential  dimension  greater  than  the 
axial  dimension  approach  those  for  a  cylinder.  Figure  5  indicates 
that,  &b  the  curvature  parameter  Z  increases,  the  critical  shear 
stresses  for  panels  having  an  axial  dimension  greater  than  the 
circumferential  dimension  deviate  more  and  more  from  the  critical 
shear  stress  for  an  infinitely  long  curved  plate.  Reference  4 
shows  that  the  reason  for  this  deviation  in  figure  5  is  that  at 
high  curvatures  the  buckling  stresses  of  these  panels,  as  veil  as 
those  of  figure  4,  approach  those  of  the  cylinder  obtained  by 
extending  the  circumferential  dimensions  of  the  panels. 
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The  effects  of  boundary  conditions  in  the  limiting  cases  of 
infinitely  long  curved  strips  (fig*  1)  and  of  complete  cylinders 
(fig.  3)  suggest  that  the  curves  of  figure  4  are  substantially 
independent  of  edge  restraint  at  large  values  of  Z  but  that  the 
curves  of  figure  5  would  be  considerably  affected  by  a  change  in 
edge  restraint. 

Long  curved  strips  tinder  combined  shear  and  direct  axial  stress . 
Reference  5  shows  that  the  theoretical  interaction  curve  fbr  a 
long  curved  strip  under  combined  shear  stress  and  direct  axial  stress 
is  approximately  parabolic  when  the  edges  are  either  simply  supported 
or  clamped,  regardless  of  the  value  of.  Z.  This  parabola  is  given 
by  the  formula 


Rs2  +  Rx  »  1 


where  Rs  and  Ex  are  the  shear -stress  and  compressive -stress 
ratios,  respectively. 

At  high  values  of  Z  curved  strips,  like  cylinders,  buckle 
at  compressive  stresses  considerably  below  the  theoretical  critical 
stresses.  In  order  to  take  this  condition  into  account,  certain 
modifications  in  the  theoretical  results  are  proposed  in  reference  5 
for  use  in  design. 

Cylinders  under  combined  shear  and  direct  axial  stress . - 
The  theoretically  determined  combinations  of  shear  stress  and  direct 
axial  stress  which  cause  a  cylinder  with  simply  supported  end 
clamped  edges  to  buckle  are  shewn  in  figure  6  (fig.  1  of  reference  6) . 
Considerable  variation  in  the  shape  of  the  interaction  curves 
occurs  for  low  values  of  Z.  For  high  values  of  Z  the  interaction 
curves  for  either  simply  supported  or  clamped  edges  are  similar  to 
the  curve  for  Z  =  30- 

Because  cylinders  actually  buckle  at  a  small  fraction  of  their 
theoretical  critical  compressive  stress,  the  theoretical  interaction 
curves  of  figure  6  cannot  be  expected  to  be  in  satisfactory  agreement 
with  experiment  when  a  very  appreciable  amount  of  compression  is 
present.  For  semiempirical  curves  and  a  check  of  available  tost 
data,  see  reference  6. 


NACA  IN  No .  1342 


17 


CONCLUDING  KEMABE3 


A  previous  investigation  showed  how  Fourier  series  type 
solutions  of  Donnell's  equation  can  he  used  to  simplify  greatly 
the  stability  analysis  of  thin  cylindrical  shells  with  simply 
supported  edges .  The  present  paper  shows  that  the  restriction  to 
simply  supported  edges  can  be  removed  by  the  introduction  of  a 
new  equation  which  is  equivalent  to  Donnell 's  equation  but  is  better 
adapted  to  solution  by  Fourier  series.  This  modified  equation  can 
be  solved  for  the  buckling  stresses  of  curved  sheet  having  either 
simply  supported  or  clamped  edges  by  established  methods  essentially 
equivalent  to  those  in  use  for  flat  sheet.  This  approach  permits 
a  simple  and  straightforward  solution  to  be  given  for  a  number  of 
problems  previously  considered  rather  formidable. 


Langley  Memorial  Aeronautical  Laboratory 

National  Advisory  Committee  for  Aeronautics 
Langley  Field,  Ya.,  March  20,  1947 
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APPENDIX 

COMPARISON  OF  RESULTS  OBTAINED  BY  USING  DONNELL'S  EQUATION 
AND  THE  MODIFIED  EQUATION  IN  THE  STABILITY  ANALYSIS 
OF  SIMPLY  SUPPORTED  CURVED  PANELS 
Solution  of  Donnell's  Equation 


Donnell's  equation  expressing  the  equilibrium  of  a  curved 
panel  under  median-surface  stresses  can  "be  written  in  general  form 
as  • 


dv8w  +  ^ 
DV  v+  r2 


+  axtV 


8x2 


(Al) 


where  x  is  the  axial  coordinate  and  y  the  circumferential 
coordinate.  Division  of  equation  (Al)  by  D  and  the  introduction 
of  the  dimensionless  stress  coefficients  kX;  k -j,  and  ks,  and 
the  curvature  parameter  Z  results  in  the  following  equation: 


_8  12Zg  sV  Jt2  4  82w  „ 

vw  +  -t-  — r  +  kY  —  V  — -  +  2kn-rV 


*2”4  ^  +  Jc„?£v4^2„0  (A2) 


"b1"  a?  & 


where 


v  n  t2t 


V  -  T't)2t 

40  T35 


and 


k  -  G 

y  ‘  Vi 


2  "  ft/1  ‘  ^ 


r 
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Equation  (A2)  can  be  represented  by 


Ql(w)  =*  0 


(A3) 


where  Qq  is  defined  as  the  operator 


d  • 


dy2 


The  equation  of  equilibrium  (equation  (A3) )  is  solved  by  using 
the  Galerkin  method  as  described  in  the  section  entitled  "Theory. " 

In  applying  this  method  the  unknown  deflection  w  is  represented  in 
terms  of  a  set  of  functions  (see  equation  (8)),  each  of  which 
satisfies  the  boundary  conditions  but  not  in  general  the  equation  of 
equilibrium..  A  suitable  set  of  functions  of  this  type,  which 
satisfies  the  boundary  conditions  for  simple  support,  is 


w 


V 


m=l 


I 


a-rrm  Bin 


MC 

a 


(A>0 


where  the  origin  is  taken  at  a  corner  of  the  plate Substituting 
in  equations  (10)  and  (ll) 


ts  sin 


mrac 

a 


&n(y)  =  sin 


and 


Q  =  Qq 


and  performing  the  integration  over  the  whole  plate  (limits  x  =  0,  a 
y  =  0,  b)  gives  the  set  of  equations 


20 


NACA  TN  No .  1342 


*  i)  * 


/ 

(i 


p<=  + 


irnpq 


(m2  -  p2)(n2  -  a2) 


0 


(A5) 


where  111=1,2,3,  n=l,2,3>  and  P  and  q  take  only  those 

value  o  for  which  m  ±  .p  and  nip  are  odd  numbers . 

Equation  (A5)  represents  an  infinite  set  of  homogeneous  linear 
equations  involving  the  unknown  deflection  coefficients  a^  ^ .  In 

order  for  the  deflection,  coefficients  to  have  values  other  than 
zero,  that  is,  in  order  for  the  panel  to  buckle,  the  determinant  of 
the  coefficients  of  the  unknown  deflection  coefficients  a.,  . 

must  vanish.  This  determinant  can  be  factored  into  two  subdeterxninants, 
one  involving  the  unknown  deflection  coefficients  a^j  for  which 

1  *  «)  is  odd  and  the  other  involving  those  coefficients  for  which 
i  i  j  is  even.  Budding  occurs,  therefore,  when  either  of  the  two 
sublets rminants  vanishes.  Only  the  buckling  criterion  involving 
the  even  subdeterminant  is  treated  here .  This  criterion  is 
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all  a13 

a22 

*31 

®33 

m=l. 

n=l 

Mi  t  0 

+^S)a 

0 

0 

m=lj 

n=3 

0  *23  - 

-5  HI)2 

0 

0 

m=2, 

n=2 

^22 

^(9+~)  += 
5vw  2 

IWEf- 

m=3. 

n=l 

0  0 

•sHt)2- 

”31 

0 

m=3, 

n=3 

r 

0  0 

+IHt)2 

0 

**33 

where 


Division,  of  each  column  of  the 
hy  the  proper 


determinant  in  equation  (A 6) 
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gives  the  simplified  equation 


The  vanishing  of  this  determinant  is  the  criterion  for  the  symmetrical 
buckling  of  the  shell.  The  same  buckling  criterion  results  from 
the  use  of  the  modified  equation,  as  is  shown  in  the  following 
section. 


Solution  of  Modified  Equation 

The  modified  equation  expressing  the  equilibrium  of  a  curved 
panel  under  median-surface  stresses  in  general  form  is 
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+  §  V4  ^  +  <rxt  |%  +  27*—^+  ffyt 


3x2 


.  a  t 

3x  3y  y  3y2 


(A8) 


Division  of  equation  (A8)  by  D  and.  simplification  of  the  result 
gives  the  following  equation: 


v4„  +  12^2  c£w  +  k  ^  +  2k  —  -5-2w  -  +  k  -  — 
hr  2t?  X  b2  3x2  s  b2  3x  3y  y  b2  3y2 


0  (A9) 


Equation  (A9)  can  he  represented  hy 


QgCw)  =  0 


(A10) 


where  Qg  is  defined  as  the  operator 


V4 


8x^ 


b2  Sx2 


£_^L+k_  £ 

h2  3x  3y  h2  Sy2 


By  use  of  the  Galerkin  method  and  hy  use  of  the  expression  for  v 
given  in  equation  (A4),  the  following  set  of  equations  analogous 
to  equations  (A5)  are  obtained 


_  _  .  laz^a1*- 

n  *) +  — 


/  2  2  a2\2 

m2  +  n2  - 


:V(m2  +  n2  J|)2 


-  ltxm2  -  hya2  ~ 


mnpq 


=  0 


+  3g£g^-  a  - 

n2h3  £_  </_  pq-  (m2  -  p2)(n2  -  q2) 

P=1  <1=1  (All) 


where  m=l,2,3;  •  ••,  n=l,2,3,  •  ••,  and  p  and  q  take  only  those 
values  such  that  m  t  p  and  n  t  q  are  odd  numbers. 
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As  In  the  case  of  the  solution  of  Donnell's  equation,  the 
stability  determinant  representing  equations  (All)  can  he  factored 
into  an  even  and  an  odd  subdeterminant.  The  even  one  is 


^l 

a13 

*22 

a3i 

a33 

m=l, 

n=l 

W11 

0 

0 

0 

m=l, 

n=3 

0 

%3 

s’ 

0 

0 

m=2, 

n=2 

4 

5 

*L2 

4 

"5 

25 

m-3 , 

n=al 

0 

0 

4 

N31 
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m=3 , 

n=3 

0 

0 

0 

n33 

2*5 


•  •  • 


=»  0  (A12) 


The  stability  determinant  (equation  (A12))  obtained  from  the 
modified  equation  is  identical  with  the  simplified  stability 
determinant  (equation  (A7))  obtained  by  use  of  Donnell's  equation. 
This  identity  holds  for  the  odd  as  well  as  the  even  determinants. 


Although  the  stability  determinants  obtained  by  use  of  the  two 
equations  are  identical  and  yield  identical  buckling  loads,  the 
determinant  in  equation  (A7)  consists  of  the  coefficients  of 


whereas  the  determinant  in  equation  (A12) 


consists  of  the  coefficients  of  a.jj.  Accordingly,  although  the 
buckling  loads  found  by  the  two  methods  are  the  same,  the  buckle 
patterns  are  different.  Of  the  two  buckle  patterns  the  one  found 
by  the  use  of  the  modified  equation  is  believed  to  be  correct. 
This  conclusion  has  been  verified  for  the  limiting  case  of  a 
flat  plate  (Z=0) . 
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TABIB  I .  -  >  WIjki  OF  PROBLEMS  TKEA325D 


Problem 

Plgara 

Reference 

Edge  Condition 

1 

E 

1 
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Singly  supported  (u-O,  Tjfo) 

Clamped  (u^O,  t»0) 

Clamped  (u»0,  Tjfc>) 

2 

2 

9 

(Leggett) 

(Simply  supported  (n«0,  rjfo) 

■1  Clamped  (ujfc>,  v-0) 
iff  lamped  (u-0,  vj to) 

(Simply  supported  (ujkt-O) 

Ifflamped  (u*v-0) 

(3D) 

3 

3 

Simply  supported 

Clamped 

■I 

m 

1 

m 

4 

4 

Singly  supported 

0 

5 

4 

Simply  supported 

1 

i 

Hot 

shewn 

5 

simply  supported 

Clamped 

1 

1 

6(a) 

6 

Simply  supported 

H 

6M 

6 

Clamped 
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Figure  2.-  Comparison  of  Leggett’s  solutions  with  present  solutions  for 
critical-shear-stress  coefficients  of  a  long  curved  strip.  (Fig.  2  of 
reference  2.) 
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.  Figure  3.-  Critical-shear -stress  coefficients  for  cylinders  in  torsion. 

(Fig.  1  of  reference  3.) 
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Figure  5.-  Critical-shear -stress  coefficients  of  simply  supported 
curved  panels  having  axial  dimension  greater  than  circumferential 


dimension,  (Dashed  curve 'estimated.)  (Fig.  2  of  reference  4.) 
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(b)  Clamped  edges. 


Figure  6.- 
stress 


Critical  combinations  of  shear-stress  and  direct-axial 
coefficients  for  cylinders.  (Fig.  1  of  reference  6.) 
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ABSTRACT: 


Modified  form  of  Donnell's  equation  for  equilibrium  of  thin  cylindrical  shells  Is  de¬ 
rived  which  is  equivalent  to  Donnell’s  equation  but  with  certain  advantages  In  physi¬ 
cal  interpretation  and  in  ease  of  solution,  particularly  in  case  of  shells  having  clamped 
edges.  Modified  equation  can  be  solved  for  buckling  stresses  of  curved  sheets  having 
either  simply  supported  or  clamped  edges,  hy  established  methods.  Question  of  im¬ 
plicit  boundary  conditions  is  also  considered. 
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